p-Wave Optical Feshbach Resonances in Yb 

Krittika Goyal, 1 - 2 ^ Iris Reichenbach, 3 and Ivan Deutsch 1 ' 2 

1 Center for Quantum Information and Control (CQuIC), University of New Mexico, Albuquerque NM 87131 

2 Department of Physics and Astronomy, University of New Mexico, Albuquerque NM 87131 
^Max Planck Institute for the Physics of Complex Systems, Nothnitzer Str. 38, D-01187 Dresden, Germany 

We study the use of an optical Feshbach resonance to modify the p-wave interaction between ultracold po- 
larized 171 Yb spin-1/2 fermions. A laser exciting two colliding atoms to the l So +^Pi channel can be detuned 
near a purely-long-range excited molecular bound state. Such an exotic molecule has an inner turning point far 
from the chemical binding region and thus three-body-recombination in the Feshbach resonance will be highly 
suppressed in contrast to that typically seen in a ground state p-wave magnetic Feshbach resonance. We calcu- 
late the excited molecular bound-state spectrum using a multichannel integration of the Schrodinger equation, 
including an external perturbation by a magnetic field. From the multichannel wave functions, we calculate the 
Feshbach resonance properties, including the modification of the elastic p-wave scattering volume and inelastic 
spontaneous scattering rate. The use of magnetic fields and selection rules for polarized light yields a highly 
controllable system. We apply this control to propose a toy model for three-color superfluidity in an optical 
lattice for spin-polarized 171 Yb, where the three colors correspond to the three spatial orbitals of the first excited 
p-band. We calculate the conditions under which tunneling and on-site interactions are comparable, at which 
point quantum critical behavior is possible. 



I. INTRODUCTION 

Alkaline-earth-like atoms are of increasing interest for 
applications in quantum control, including optical atom 
clocks IT), quantum computing 1 2 -5 1, and simulations of con- 
densed matter systems (6). Experimental advances are pro- 
ceeding at a steady pace with demonstrations of a variety of 
important milestones, including clocks that now surpass the 
cesium standard [7|, Bose-Einstein condensation in isotopes 
of ytterbium [8], calcium [9|, and strontium ifTUl [Til . Fermi 
degenerate gases lfT2l . and the superfluid-to-mott-insulator 
quantum phase transition ifTSIl . 

Another important ingredient in the quantum-control tool- 
box is the ability to control the interatomic interactions. Fes- 
hbach resonances have played an essential role in such ma- 
nipulation of alkali-metal degenerate gases, allowing for the 
observation of the BEC-BCS crossover [14J. Whereas in al- 
kali gases Feshbach resonances can be induced via magnetic 
fields that couple different channels in the electronic ground 
state, in alkaline-earth-like atoms this is not possible because 
of the lack of hyperfine structure in the ground 1 So state. An 
alternative is to employ an optical Feshbach resonance (OFR) 
by laser-coupling two scattering ground-state atoms to a meta- 
stable bound molecule in an excited-state potential [ 15 1. Alka- 
line earths are particularly well suited to OFRs due to the exis- 
tence of narrow intercombination lines of the kind studied for 
optical clocks, 'So — > 3 Pj [ 16 1. Photoassociation spectroscopy 
has been used to measure narrow molecular resonances in the 
'So + channel lfT7l[T8l . an important first step toward im- 
plementation of OFRs. 

In previous work we studied the use of OFRs to manip- 
ulate nuclear spin coherence in fermionic, spin- 1/2, 171 Yb 
using OFRs associated with s-wave collisions |[T9l . In the 
work presented here, we extend our study to p-wave OFRs 
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of this species. The ability to manipulate p-wave collisions 
could open the door to studies of nonconventional superflu- 
idity and other exotic quantum phases of matter [20|. Prior 
observations of p-wave magnetic Feshbach resonances in al- 
kalis proved to be too lossy for quantum coherent control ETI . 
Inelastic collisions are believed to be enhanced in these reso- 
nances because the p-wave scattering states are well localized 
behind the centrifugal barrier ED . They thus have a very large 
Franck-Condon overlap with more tightly bound molecules 
below the Feshbach threshold, which leads to exothermic tran- 
sitions. The use of OFRs can potentially mitigate this effect. 
In particular, purely-long-range (PLR) molecular states ex- 
isting in excited state potentials can be coupled optically to 
ground-state p-wave channels [18|. These PLR states, arising 
from avoided crossings in the excited state hyperfine structure, 
have inner turning points at ~ 50ao, and are thus well sepa- 
rated from the chemical binding region. Inelastic collisions to 
bound-ground molecules via excitation to PLR states should 
be highly suppressed. In this case, heating due to spontaneous 
emission will be the dominant source of inelastic collisions, 
but this too can be suppressed through off-resonance excita- 
tion. 

In addition to suppressing inelastic recombination, OFRs 
offer opportunities for quantum control beyond what is possi- 
ble with magnetic Feshbach resonances. For example, in p- 
wave collisions the projection of rotational angular momen- 
tum along a given axis is a new degree of freedom that can 
affect the symmetry of the order parameter in p-wave super- 
fluidity ll22l . In the presence of a bias magnetic field and for 
appropriate choices of laser polarization, we can address these 
degrees of freedom and control the scattering length associ- 
ated with different projection quantum numbers. If an optical 
lattice trapping potential is added, a variety of rich phenom- 
ena can be explored with such control. For example, the three 
projections of angular momentum translate into three orbitals 
of a p-band in the first excited vibrational state of an optical 
lattice E3l l24l . With control of p-wave collisions of spin- 
polarized fermions, one can obtain a Hubbard model similar 
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to the one that gives rise to 3-color superfluidity and trionic 
phases and is an important model of QCD E51 . 

In this article we study the use of an OFR to control p- 
wave collisions in 171 Yb by exciting near photoassociation 
resonances of the 1 Sq + 3 Pi channel. After reviewing the sys- 
tem and the formalism for calculating the optically controlled 
scattering properties, we calculate the energy spectrum and 
scattering lengths including the presence of a magnetic field, 
which allows for polarization-dependent control of the inter- 
action. We apply this to a toy model of 3-color superfluidity 
to give a benchmark of the performance of the p-wave OFR 
and summarize our results. 
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II. p-WAVE PHOTOASSOCIATION RESONANCES 

We consider spin-polarized 171 Yb, with nuclear spin ; = 
1 /2, for which s-wave collisions are forbidden and p-waves 
dominate at low temperature. The essential formalism for de- 
scribing the system, in the absence of an external magnetic 
field, was given in fl9l . We review the salient points here. 
The two-atom states in each of the collision channels are gov- 
erned by an effective potential of the form 



V ef f = 
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where Vbo is the Bom-Oppenheimer potential in Hund's case- 
(c), Vhf is the hyperfine interaction, and V mag is the interaction 
with external magnetic fields. Here and throughout we set 
h = 1 and we use atomic units. For the ground So + 1 Sq col- 
lision there is only one channel, the nuclear spin triplet state 
/ = 1 , mi = 1 . There is no hyperfine interaction and we neglect 
the very small magnetic interaction with the nuclear magne- 
ton. As we are interested only in the near-threshold scat- 
tering states of this channel, the ground Born-Oppenheimer 
potentials can be approximated in a modified Leonard-Jones 
form 1261 . 
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where C { 6 8) = 1931. 7a.u., = 1.93 X 10 5 a.u., and C$ = 
1.03409 x 10 9 a.u. l26l . Since we are considering p-wave 
scattering, the rotational angular momentum is R = 1. The 
system is not prepared in a state with a fixed projection of R, 
and thus the atoms can scatter with any allowed value of mr = 
— 1,0, 1 relative to a space-fixed quantization axis, defined by 
the magnetic field. We obtain the scattering wave functions 
corresponding to the above potential numerically, using the 
Numerov method for integration 6 271 l28l . 

In the excited 1 Sq + ^Pi channel, the description is more 
complicated. The electronic Born-Oppenheimer potentials are 
taken in the Hund's case-(c), 
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with parameters determined by fits to experiments as Ci = 
2810a.u., C$ = 1.862 x 10 8 a.u. and Cf =1 = -Cf =0 /2 = 



Figure 1 . Adiabatic potentials for the four channels with T = 3 that 
asymptote to the 1 So + {fl = 3/2) channel. Since Mj takes seven 
values, each channel is seven fold degenerate. 



0. 09695a. u. for the 1„ and 0„ states respectively. The Hund's 
case-(c) variables, however, are not good quantum numbers in 
the region of interest. Coriolis forces mix nuclear rotation and 
electronic angular momentum and hyperfine interaction mixes 
this with nuclear spin [29]. As such, the only good quantum 
numbers are the total angular momentum and its projection 
which we denote T,Mj; parity is fixed here to be -1 for the 
p-wave collisions. 

We are interested in the molecular bound states, or photoas- 
sociation resonances of these electronic potentials. Dipole se- 
lection rules break the resonances into two parity classes - 
those accessible from P-even or /?-odd ground states |[T8l[T9l . 
Of particular interest are the PLR states arising from avoided 
crossings due to hyperfine mixing. Figure [T] shows the adi- 
abatic potentials with T = 3 that asymptote to the 1 Sq + 
3 P\(f2 = 3/2) channel, where fi is the hyperfine quantum 
number of the excited state atom. There exists one poten- 
tial with its minimum at ~ 75ao and a depth of 0.68 GHz. 
This shallow potential nonetheless supports bound states that 
are well resolved and can be used for p-wave OFRs with sup- 
pressed three-body recombination. 

To determine the photoassociation resonances, we employ 
a multichannel integration of the Schrodinger equation as dis- 
cussed in |fl9l . We consider first the case of no external mag- 
netic fields. The effective potential operator in the 1 5'o + 3 Pi 
channel, Eq. ([lj, is written as a matrix expanded in the ex- 
tended Hund's case-(e) basis \e(T,M T )) = \f 2 ,F,R,T,M T ), 
where F = fj + f 2 and T = F + R (39). Here /i = 1/2 is 
the spin of the ground-state atom, and f% = 3/2,1/2 is the 
hyperfine spin of the excited-state atom. In the ground state 
I = F = 1, Mf = 1, R = 1, and the total angular momen- 
tum takes the possible values T g = 0, 1,2. By dipole selection 
rules, in the excited channels the allowed values are therefore 
T e = 0, 1,2,3. The effective excited potential matrix thus has 
19 channels each of which are 2T e + 1 fold degenerate, re- 
sulting in a total of 89 channels. We denote the multichannel 
excited bound states as (neglecting the subscript e), 

\n,T,M T }= £ y B ^ w (r)|e(r,Jlfr)). (4) 

e(r,M r ) 
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Figure 2. Spinor components of the multichannel wave function of 
the PLR bound molecular state at -355 MHz and with T = 3. Each 
curve corresponds to a component associated with one of the six ba- 
sis states, \z(T,Mt)}, that contribute to this state. Since Mj can as- 
sume seven distinct values, each wave function is 7-fold degenerate. 



In the binding energy range of —1022 MHz to —3 MHz, the 
system supports 2 bound states with T = 0, 26 bound states 
with T = 1, 15 bound states with T = 2, and 23 bound states 
with 7 = 3. 

Of particular interest are the PLR states, denoted in Table 
[l|a). Figure [2] shows an example of a multichannel spinor 
wave function of the PLR bound molecular state at —355 MHz 
with T = 3. Each spinor component corresponds to one of the 
six different \e(T,Mx)) channels, each of which are 7-fold 
degenerate. Most of the amplitude of the wave function is 
supported between 50ao and 150ao- As such, the inner turn- 
ing point is well removed from the chemical binding region 



and the outer turning point is sufficiently far out to allow for 
a large Franck-Condon factor in optical excitation. These fea- 
tures are advantageous for application to OFRs. 



A. In an external magnetic field 

We now consider the effect of an external magnetic field 
to allow for additional control on the system. With the B- 
field defining the quantization axis and in the linear Zeeman 
regime, the perturbing potential is 



Vmag= £ 8f 2 /JBBm f2 \f2,m f2 )(f2,m f2 

f2,>"f 2 
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where gf. is the Lande g-factor of the atomic hyperfine level. 
This Hamiltonian breaks the rotational symmetry and gener- 
ally couples an infinite hierarchy of states with different total 
angular momenta T. For the relatively weak magnetic fields 
that we consider here, we can employ perturbation theory. We 
break the degeneracy of the states within a r-manifold and 
mix states with the same Mj when the Zeeman shift is on 
the order of the vibrational spacing. For these weak magnetic 
fields, the value of T at zero magnetic field still dominates and 
this will be used to label the states. This is particularly true for 
the PLR states, where T remains approximately a good quan- 
tum number for all fields we use in our calculation. 

To obtain the eigenenergies and eigenfunctions in the mag- 
netic field, we diagonalize V mag expressed as a matrix in the 
basis of the bound states \n,T,Mr) within the energy range 
given in the discussion following Eq. (4). The matrix elements 
are given by 



(»,r,M r |y mag |«',r',Mi 
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The term (z(T ,M T )\V mag \z' (T ,M T )) characterizes the cou- 
pling of the spin degrees of freedom and the Franck-Condon 
overlap, J\\f* £{TMT) (r)\\f n ,^ {rMT) (r)dr, is the coupling of 
the radial wave functions. 



A part of the eigenspectrum, between -427 MHz and -273 
MHz, is shown in Fig. [3] The PLR state of interest, with bind- 
ing of 355 MHz and T — 3, exhibits an approximately linear 
Zeeman splitting of its 7 magnetic sublevels over a range of 
80 Gauss, as shown in the inset. Figure [3] also shows that the 
PLR state, with binding of 383 MHz and T = 1, also has an 
approximately linear Zeeman splitting of its 3 magnetic sub- 
levels over the 80 Gauss range, while the remaining two states 
(-279 MHz and -416 MHz), which are not PLR, show nonlin- 
ear Zeeman shifts over that range of perturbation. 



III. THE p-WAVE OFR 



To calculate the effect of the OFR on the p-wave scattering 
volume we turn to the theory of Bohn and Julienne ll30l . In 
that formalism the laser field is chosen detuned close to, but 
off -resonance from, a given photoassociation resonance. Only 
one bound state in a closed channel is assumed to contribute 
to the modification of the scattering volume. In practice, the 
laser field can couple to multiple excited bound states, and 
in the far-off-resonance limit, all will contribute. How such 
multiple resonances interfere and affect the OFR is a subject 
of continued research. Here, we will choose parameters for 
which one PLR bound state dominates and calculate its con- 
tribution to the OFR in both elastic and inelastic terms. 

For a single bound state, the effect of the OFR on the S- 
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Figure 3. Eigenspectrum of states in Table[IJb) (energy range: -427 
MHz to -273 MHz) as a function of an applied magnetic field, B. 
For B = the eigenenergies and the quantum number T correspond 
to the values in the left most column of Table |l{b). We calculate 
the OFR associated with tuning near the T = 3 PLR state, bound 
by —355 MHz, which splits into 7 magnetic sublevels in a linear 
Zeeman regime of the 80 Gauss plotted here 



matrix in the incoming So + Sq channel is 



2 ,, n 2A-/(r-y) 
2A + i(r + y)" 



(7) 



y is the molecular natural linewidth, A is the detuning of the 
laser from the bound molecular state (including the light-shift 
of that level), rjo is the background phase shift, and the stimu- 
lated linewidth is, 
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Here / is the laser intensity, I sat =0.13 mW/cm 2 is the atomic 
saturation intensity for 3 Pi, and Ya/271 = 182 kHz is the 
atomic linewidth. The Franck-Condon factor, with rotational 
corrections, is 
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expressed here as the ratio of the free-to-bound transi- 
tion molecular dipole moment for laser polarization El to 
the atomic dipole momentum d\ = (3c 3 y4)/(4co 3 ). In the 
Wigner-threshold regime |Vg(^r)| 2 05 where k r is the wave 
vector of the relative coordinate momentum at the scattering 



energy, and thus we define 1/ opt as the "optical volume" in 
analogy with the "optical length" for s-wave OFRs, 
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This is the parameter that defines the strength of the /j-wave 
OFR. 

Selection rules dictate the allowed transitions that are ac- 
cessible for an OFR. For atoms with spin-polarized nuclei 
scattering on the ground 1 Sq + 1 So potential, 1=1, mi = 1, 
T g = 0, 1,2, and Mj g = yiir + 1, where mg = — 1,0, 1 are the 
projections of the partial-wave angular momentum on the 
quantization axis. We can thus optically connect to excited 
Sq + 3 P\ bound molecules with T e = 0,1,2,3, and Mj e = 
Mt. +q = mn+l+q, where q denotes the projection of pho- 
ton helicity (n, a± ). Table[ljb) shows the values of % pt for the 
coupling of the partial- wave projections to four of the ex- 
cited bound molecular states using a+ polarized light at B = 
Gauss and B = 30 Gauss. A non-zero B-field leads to a mix- 
ing between the different eigenstates with the same Mj e , lead- 
ing to a change in the Condon overlap of this state with the 
scattering wave function of the ground potential. For the PLR 
states we see that the T^, pt is fairly constant, while for the other 
states the V 0?t is significantly changed by the magnetic field. 
This is because the poor overlap of the PLR states with their 
neighboring non-PLR states suppresses mixing. 

With the scattering matrix in hand, the p-wave scattering 
volume is defined as a 3 = —K/k^, where the Zf-matrix ele- 
ment is given by 13D 
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excluding the background phase shift. The real and imaginary 
parts of the p-wave scattering volume are then 
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where the additional background contribution was added. We 
obtain the background phase shift r|o by numerical integration 
of the p-wave scattering state in the Wigner threshold regime 
and fit to the asymptotic wave function. We find the back- 
ground scattering volume to be a\ g = —406446 a.u. 

The real and imaginary parts of the scattering volume gov- 
ern the strengths of the elastic and inelastic collisions respec- 
tively. In principle, one can increase the ratio of good to bad 
collisions solely by increasing the detuning. In practice, this 
is limited by the available intensity that is required to ensure 
a sufficiently strong interaction. Moreover, our model is re- 
stricted to an OFR via a single molecular bound state, and for 
self-consistency, we require a sufficiently small detuning so 
that only one photoassociation resonance dominates the pro- 
cess. For these reasons, we must choose a state in a suffi- 
ciently sparse region of the density of states so that when the 
laser is detuned closest to this state, even for detunings large 
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Table I. p-wave optical volumes (T^opt) for the coupling of all possible initial states to four of the bound molecular states of the excited potential 
(the energies of these states and their T value are shown in the first column ) with (a) different polarizations, (b) polarization q = 1. Blank 
entries indicate that the particular combination of initial state, polarization and final state is forbidden. A * indicates that the particular state is 
a PLR state. The PLR state bound at —355 MHz, denoted in bold face, is used for the OFR calculation presented here. 



enough to avoid spontaneous scattering, the single resonance 
model is valid. 

We thus seek a PLR state that we can address with high res- 
olution and with a sufficient optical volume to induce a strong 
OFR. Firstly, as we are considering spin-polarized fermions, 
we can ignore the nearby spectrum of s-wave photoassocia- 
tion resonances and concentrate only on the bound states con- 
nected to /?-waves. Secondly, by employing dipole selection 
rules, we can reduce the number of allowed transitions and 
reduce the density of states. Using a magnetic field and po- 
larized light, the interaction strength for scattering in states 
of the ground potential with a particular m R value can be se- 
lectively enhanced while suppressing the interaction strength 
for scattering in states with other m R values. For example, 
the state with m R = 1, corresponding to the stretched state, 
r ? = 2,Mj g = 2, couples with o + polarized light only to a 
T e = 3,Mj e = 3 molecular bound state. Transition to states 
with other values of T e are forbidden. Of course, the ground 
state can not be prepared in a state with a given m R , but in 
the presence of a magnetic field, differences in detuning and 
optical volumes can suppress other transitions. 

Given these observations, the PLR at —355 MHz is promis- 
ing for application to a p-wave OFR. This is a T e = 3 state 
which connects only to a T g = 2 ground-state. In the presence 
of a magnetic field and with o+ polarized light, we can ad- 
dress the ground Mj = 2 — > Mj e = 3 transition and make this 
the dominant resonance (see Fig. |4}. The neighboring bound 
states are T e = \ (see Fig. 3) and inaccessible with this polar- 
ization from the T g = 2,Mj g = 2 ground state. In addition, the 
T^pt for m R = 1 is substantially larger than V 0?x for the other 
sublevels indicating that the OFR is strongest for the Mj g = 2 
state. This leads to a further enhancement of the interaction 
strength for the Mr = 2 — > Mj e = 3 transition. 

Figure |4] shows a possible configuration for inducing the 



OFR. In a 30 Gauss magnetic field, and detuning A = — 3 MHz 
below the resonance at —355 MHz, we dominantly couple the 
T g = 2.M Tg = 2 -> T e = 3,M Te = 3 transition. Using Eq. (fl2|) 



we calculate the real-part of the scattering volume arising 
from the OFR to be 3i(a 3 p ) = -1.44 x lO^W/cm 2 )- 1 . The 
imaginary part is reduced by the factor j/2A = 0.057. The 
effect of this spontaneous emission will depend on the appli- 
cation at hand. Coupling to other transitions, Mj g = 0, 1 — > 

M Te = 1,2 are reduced to ^.{a 3 p ) = -4.36 x 10 4 (W/cm 2 )- 1 

and 9t(4) = -1.40 x lO^W/cm 2 )- 1 respectively. In addi- 
tion, off -resonant coupling of T g = 2,Mj g = to neighboring 
T e = \ .Mj e = 1 is highly suppressed at this detuning. 



IV. THE THREE-COLOR FERMI-HUBBARD MODEL 

It will be extremely challenging to observe /?-wave super- 
fluidity in a dilute gas, even with the use of an OFR, given the 
ultra-low temperatures required. Nonetheless, the ability to 
control p-wave interactions can potentially lead to a rich va- 
riety of many-body phenomena, particularly if an optical lat- 
tice confining potential is included. We propose here how the 
combination of such tools can be used to explore a toy model 
of fermionic color superfluidity with three colors. Such mod- 
els have been considered before 1251 where the internal de- 
grees of freedom served as the three "colors". For the case of 
fermions, this is not a natural realization since the number of 
internal states will always be an even number. An alternative 
is to employ the external degrees of freedom associated with 
the three spatial orbitals of the first excited "p-band" of an 
optical lattice. Such colors have been considered for bosons, 
mediated by s-wave interactions. We consider here a model 
for spin polarized fermions, mediated by p-wave interactions. 

Following [24], the multicolor field operator for spinless 
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(i.e. polarized) fermions in the first excited /?-band is written 
in the Wannier basis as 



\|/(x) =^c,- a (j) a (x-R ; ), 



(14) 



where <^ a (x) is a p-orbital with a = x,y,z, and c, a is the 
fermionic annihilation operator for that orbital at the i th lat- 
tice site. We consider lattices of sufficient depth Vq that the 
tight-binding approximation is valid. We restrict the dynam- 
ics to a single p-band, which can be metastable, as seen in 
recent experiments where bosons remained in the first excited 
band of an optical lattice for about a hundred times the tun- 
neling time scale ll32l[33l . We expect a similar metastability 
for fermions. In addition, we assume sufficiently deep lattices 
such that the tunneling coefficient for a particle in the state a 
in negligible along the direction a' for a' 7^ a. Moreover, we 
take the wells to be spherically symmetric. The Hamiltonian 
then takes the Fermi-Hubbard form for the three colors in a 
single band, 

c la c j-a+ H c i,a' c lf,' Ci V c i.?> V ai$',a$> 
(i-j)aA i,m«,ap,a'|3' 

(15) 

where J is the tunneling coefficient along any direction a, 
(i,j)a indicates that i and j are nearest neighbors along a, 
and V^"p, a p is the interaction matrix element for two atoms at 
the same site starting in orbitals a, p and scattering to oc',p' 
via p-wave collisions of symmetry niR. The coupling matrix 
is 



H= J £ 



v mR 



<&( Xl )4£( X2 )V"*(xi -x 2 )4>a(xi)<|>p(x 2 ), (16) 



where V™ R 



(xi — x 2 ) is the two-body interaction potential for 
p-wave scattering. This can be treated through a pseudopo- 
tential on a delta-shell l34l 



33i(al) 



6(r-j) 3 3 



(17) 



In order to calculate the interaction matrix we transform the 
Wannier states from the Cartesian orbitals to spherically sym- 
metric 3D harmonic oscillator orbitals, and to center-of-mass 
and relative coordinates of the two particles, specified by the 
projections of angular momentum, Mr and ihr, respectively. 
The matrix then takes the form 



V mR 



■ Y,(u'V\m r M R )U m *(m r M R \aV), (18) 



Mr 



where (m r Mn |ap) is the angular part of the change-of-basis 
matrix, and U mR is the interaction strength coming from the 
radial integral of the interaction potential expressed in the rel- 
ative coordinate, proportional to the real-part of the p-wave 
scattering volume. 

Like the model studied in ll25ll . the Fermi-Hubbard Hamil- 
tonian Eq. ( 15 1 has three colors, but differs in two important 



ways. Firstly, it allows for anisotropic interactions as consid- 
ered in [35 1 . In addition, we allow for couplings between dif- 
ferent incoming and outgoing orbitals, a 7^ a', p 7^ p', as stud- 
ied for bosons, ll24l . Most importantly, unlike any model pre- 
viously considered, the control provided by the OFR allows 



Mj = 1 



, , Mj 
3 ^--2.18x 10 6 a0 3 — 



-6.97 x 10 5 aO 3 

»(4 




3 MHz 



mr — — 1 niR = niR = 1 
Mr = Mr. = 1 Mj = 2 



Figure 4. OFR using a + polarized light to couple the scattering 
state of the ground potential with the three different projections of 
p-wave angular momentum, mr, to the excited PLR bound state with 
the different total projection Mj e . The figure shows only those states 
permitted by selection rules. Denoted are OFR values of the real part 
of the scattering volume, 5K(a 3 ), for each of the three transitions, for 
an intensity of 50 W/cm 2 . 



for the possibility to manipulate the strength of interactions 
in a manner that depends on the fermionic colors. We expect 
such control could be used to explore a variety of phenomena 
such as the trionic phase and color superfluids discussed in 
Ref. Il25l . We leave the details of the many-body analysis for 
future work. 

To evaluate the potential for this system to lead to quantum 
critical behavior, we give here a rough back-of-the-envelope 
estimate. We expect that interesting many-body physics will 
be accessible when the ratio between the kinetic and interac- 
tion energies in the system is of order one [36], i. e., U mR > 67, 
for some given choice of rag. Choosing the lattice depth along 
any direction to be Vo = 18£V, where E r is the recoil energy, 
we find J = 0.16£V in the first excited band. For these param- 
eters, it follows that phase transitions occur near 



|^(a 3 )|>7x 10 6 a.u. 



(19) 



Typically, the p-wave scattering volume arising from the 
background phase shift is very small and the model in Eq. ( p"5j ) 
does not result in quantum phase transitions. However, us- 
ing an OFR, 5R(a 3 ) can be tuned to larger values. Moreover, 
through selection rules we can control specific m^-couplings 
that correlate with interactions of specific colors. Figure [4] 
outlines one possible scheme. The different Mj g levels are 
coupled to specific Mj e levels in the excited PLR state, shown 
in Fig. [3] using a + polarized light. The Zeeman splitting be- 
tween the different Mj e levels of the excited state is approxi- 
mately 0.89 MHz at B = 30 Gauss. The figure indicates the 
values of 9i(a 3 ), calculated for the couplings between the dif- 
ferent states for a laser of intensity / = 50W/cm 2 and a detun- 
ing A = 3 MHz below the Mj e = 3 state. For this magnetic 
field, laser intensity, polarization, and detuning the atoms, 
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scattering in the m# = 1 state will experience a p-wave scatter- 
ing volume of 5R(flp) = —7.24 x 10 6 a.u., satisfying the crite- 
rion in Eq. (19) . With such control, we expect one can observe 
novel quantum critical behavior in the fermionic superfluid. 

V. SUMMARY AND OUTLOOK 

We have studied a highly controllable system of spin- 
polarized 17 *Yb atoms undergoing p-wave collisions as mod- 
ified by an optical Feshbach resonance (OFR). By tuning near 
an electronically excited purely-long-range (PLR) bound state 
in the 1 So + 3 A channel, we expect to suppress three-body re- 
combination losses that typify magnetically induced /j-wave 
Feshbach resonances in the ground electronic manifold. We 
used a multichannel integration of the Schrodinger equation to 
determine the photoassociation resonances and the eigenfunc- 
tions including perturbing magnetic fields. With these, we cal- 
culated the real and imaginary part of the "scattering volume" 
associated with the p-wave scattering phase shift and loss rate 
for choices of magnetic fields and OFR polarized laser fields. 

Because the 3 Pi state has a relatively large linewidth as 
compared to the other intercombination lines, the demands 
on precision control of laser detuning are moderate. On the 
other hand, this larger linewidth implies a limitation on the 
strength of the real-part of the Feshbach resonance before in- 
elastic scattering can no longer be neglected. For these rea- 



sons we expect that even with an OFR, one will not be able to 
achieve p-wave superfluidity, or a BEC-BCS crossover, anal- 
ogous to that seen for s-wave pairing. Nonetheless, the degree 
of control afforded by the OFR could open the door to explo- 
rations of novel quantum critical behavior in the many-body 
system. 

We began such an exploration, considering a new model of 
three-color fermonic superfluidity. Here the three colors cor- 
respond to the three spatial orbitals of spinless (i.e. polarized) 
fermions in the first excited p-band of an optical lattice. Based 
on this toy model, we calculated the parameters of a Hubbard 
model including nearest neighbor hopping and on-site inter- 
action between two fermions in different orbitals via p-wave 
collisions. Through careful choice of magnetic field, laser po- 
larization, and detuning, we find conditions under which tun- 
neling and interaction energy scales are comparable. For such 
operating conditions, we expect quantum phase transitions are 
possible. A full many-body exploration of the phase diagram 
is left for future analysis. 

We thank Maciej Lewenstein, Pietro Massignan, and 
Philipp Hauke for helpful discussions, particularly about the 
application of our model to the Fermi-Hubbard Hamiltonian. 
KG and IHD acknowledge support from the Office of Naval 
Research Grant No. N00014-03-1-0508 and the Center for 
Quantum Information and Control (CQuIC) via the National 
Science Foundation Grant PHY-0969997. 
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